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I ourier  analysis  techniques,  provides  output  in  the  form  of  Doppler  frequency 
vs  arrival-angle  maps  as  well  as  coherence  maps.  These  maps  are  used  to 
analyze  effects  of  the  irregular  high-latitude  ionosphere,  such  as  Doppler  fre- 
quency and  arri\ al-angle  spreading  of  the  signal,  and  the  temporal  and  spacial 
coherence.  Various  antenna  configurations  and  spacings,  together  with  adap- 
tive processing  techniques,  permit  optimization  of  he  resolution  and  signal/ 
noise  ratio  of  the  various  parameters  being  analyzed.  Use  is  made  of  bistatic 
oblique  ionograms  and  a special  polar-model  computer  prediction  program  to 
compliment  the  analysis. 

This  report  presents  selected  data  and  the  initial  results  of  the  effects  on 
forward  propagation  using  the  DA  AS  \1  system  in  conjunction  with  a moving  ail — 
craft  to  provide  signal  Doppler  and  arrival-angle  information.  A KC-135 
specially  equipped  for  ionospheric  research  made  six  8-hour  flights  of  up  to 
3000  km  into  the  Arctic  from  Goose  Bay. 

Examples  of  data  selected  to  illustrate  particular  propagation  effects  are 
presented,  such  as  Doppler  frequency  spreading,  arrival-angle  spreading, 
multiple  mode  propagation  causing  deviation  from  predicted  values  of  both 
arrival-angles  and  Doppler  frequencies,  etc. 

Two  of  the  most  serious  problems  indicated  by  the  Goose  Bav  DA  ASM 
experiment,  in  terms  of  the  operation  of  OTH  radar  systems,  appear  to  be 
the  effects  of  TID’s  and  tilted  ionospheres.  Investigation  of  the  problem  is 
continuing  and  further  results  of  the  analysis  will  be  presented  in  a forthcom- 
ing report. 
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DAASM  Project  - High  Latitude 
Aircraft  HF  Propagation  Experiment 


1.  INTROUl  C.TION 


A series  of  papers'  described  n method  of  vertical  ionospheric  "drift"  mea- 
surements and  analysis  using  complex  amplitude  measurements  at  spaced  re- 
ceivers. This  method  o'  analysis  first  Fourier  analyzes  the  received  signal  at 
sach  antenna  into  its  Doppler  frequency  components,  end  then  for  each  Fourier 
component,  determines  the  angle  of  arrival  using  the  spatial  array  of  antennas. 
This  technique  has  been  adapted  into  an  oblique  sounding  system  using  a linear 
array  of  antennas  and  designed  to  operate  either  on  hackscatter  signals  or  signals 
propagated  from  a similarly  equipped  aircraft  in  ; forward-propagation  mode.  The 
DAASM  system  (DAASM  is  ar,  acronym  for  Doppler/Arrival  Angle  Spectral  Mea- 
surements) has  operated  at  Goose  Hay,  Labrador  (5.^20'N.  60°20'\V)  for  just  one 
and  one-half  years  beginning  full  operation  in  January  1074. 

(Received  for  publication  10  May  10751 

1.  Pfister,  W.  , et  al  (1068-75)  Pulse  Sounding  with  Closely  Spaced  Receivers  as 
as  a Tool  for  Measuring  Atmospheric  Motions  and  Fine  Structure  in  the 

Ionosphere,  ’ ___ 

VoTTTT  Environmental  Research  Paper  No.  295,  Dec  1-68 
Vol.II,  Environmental  Research  Paper  No.  T95,  Dec  1968 

Vol.III.  Environmental  Research  Paper  No.  3 17 , Mar  10^0 

Vol.IV,  Environmental  Research  Paper  No.  329,  Aug  1970 
Vol.  V,  Environmental  Research  Paper  No.  468,  F eb  1974 

Vol.VI,  Environmental  Research  Paper  No.  47  0,  Mar  1974 

Vol.  VII,  Environment -1  Research  Paper  No.  506,  Apr  1975 
Vol.  V HI,  Environmental  Research  Paper  No.  507,  Apr  1975 
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The  purpose  of  the  l)A ASM  program  is  to  investigate  the  effects  of  the  auroral 
and  polar  ionosphere  on  a backscatter  sounding  system  which  at  the  same  time  can 
receive  a forward-propagated  signal  from  an  airborne  transmitter.  The  aircraft 
used  for  these  experiments  was  a KC-135.  specially  equipped  for  ionospheric  ex- 
periments. The  questions  to  he  answered  concerned  the  character  of  the  received 
signals  and  are  focussed  on  three  main  areas: 

<al  Angle  of  Arrival  - The  degree  to  which  the  signals  from  the  aircraft  is  ' 
deviated  from  the  known  bearing  to  the  aircraft  and  the  amount  of  angular  spread- 
ing associated  with  the  particular  path. 

(b)  Doppler  Frequency  - The  extent  to  which  propagation  jn  the  auroral  re- 
gions deviates  the  expected  Doppler  frequency  associate,!  with  the  aircraft  motion 
and  the  amount  of  !r-,,ucnrv  -neading  that  occurs  ..n  the  Mgnal. 

(cl  C oherence  - The  temporal  and  spatial  coherence  of  these  aircraft  signals 
after  propagation  through  the  irregular  auroral  medium. 

The  basic  concept  of  the  DA  ASM  system  is  to  record  on  digital  magnetic  tape 
the  signal  output  of  each  antenna  in  the  receiving  array.  This  recorded  data  repre- 
sents a two-dimensional  description  jf  the  received  wave  field,  that  is,  time  and  a 
one-dimensional  space  (linear  antenna  arm"1.  In  the  most  simple  approach,  a 
two-dimensional  Fourier  transform  of  this  data  set  results  in  a Doppler  frequency 
vs  arrival  -angle  description  of  the  received  signal.  These  transforms  are  accom- 
plished at  the  AFC-HL  computer  facility  after  the  tapes  are  returned  from  the 
Goose  Hay  site. 

I he  sampled  data  is  limited  both  in  time  and  space  because  of  the  finite  dura- 
tion of  the  sample  (32  sec)  and  the  limited  extent  of  the  antenna  array  (DIO  m 
maximum).  The  resulting  Doppler  frequency  and  arrival  angle  description  of 
the  signal  is  thereby  limited  in  resolution.  This  is  particularly  true  of  the  spatial 
dimension,  and  advanced  processing  techniques  having  been  developed  to  improve 
the  capability  of  the  system  in  this  area. 

This  report  addresses  the  six  aircraft  missions  carried  out  during  11)74.  and 
an  attempt  is  made  to  relate  the  observed  character  of  the  received  signals  to  the 

propagation  mode  and  other  geophysical  factors  which  may  affect  the  auroral  iono- 
sphere. 

2.  KXPKRIMKNT 


2.1  I;  A ASM 

The  DA  ASM  system  as  configured  for  the  Goose  Bay  experiments  measures 
the  received  amplitude  and  phase  of  an  arriving  radio  signal  as  a function  of  both 


time  and  space,  in  order  to  accomplish  this,  an  array  of  small  receiving  antennas 
was  erected  in  a linear  configuration  broadside  to  the  required  receiving  direction. 
In  this  case  the  boresight  of  the  array  at  (loose  Bay,  that  is.  the  direction  normal 
to  the  length  of  the  array,  was  magnetic  North,  328°  T.  In  order  to  provide  maxi- 
mum flexibility  in  terms  of  spatial  sampling,  a complex  antenna  arrangement  was 
2 3 

developed  using  a 20-element  array.  The  DAASM  electronic  system  was  de- 
signed to  use  subsets  of  this  array,  consisting  of  either  12  or  6 elements.  Mach 

a 

subset  has  either  a uniform  spacing  or  a Minimum  Redundancy  Array  (MRA)/' 

The  MBA  provides  all  spacings  that  are  multiples  of  a basic  spacing  up  tc  some 
maximum  spacing  with  the  minimum  number  of  antennas.  In  the  DAASM  system 
the  several  antenna  sub-array  configurations  available  are  briefly  described  in 
the  following  table: 


Table  1.  DAASM  Arrays  (Large) 


No.  of 
Elements 

Min.  Spacing 

Max.  Spacing 

Remarks 

Configu  ration 
Code 

12 

1 0 m 

110  m 

Uniform 

5/0 

1 1 

1 0 m 

4 50  m 

MRA 

5/4 



The  configuration  code  above  is  designated  by  the  symbols  T/A,  where  T 
selects  the  frequency  mode  and  sampling  sequence  and  A selects  the  antenna  se- 
quence. As  will  be  described  later,  the  DAASM  system,  was  designed  to  operate 
with  pulsed  transmissions  at  a 400  liz  repetition  frequency.  The  experiments 
described  in  this  report  all  involve  measurements  carried  out  in  conjunction  with 
the  transmitter  on-  board  the  KC-135  aircraft  as  a signal  source  flying  to  the 
north  of  Goose  Bay.  Because  of  transmitter  average  power  limitations  on  the 
aircraft,  pulse  repetition  frequency  was  limited  to  200  Hz.  This  necessitated  a 
change  in  the  Goose  Bay  antenna  system.  Five  additional  sub-arrays,  each  with 
six  elements,  were  devised.  Table  2 describes  these  sub-arrays.  All  of  the 
arrays  are  included  in  the  20  antennas  and  are  selectable  by  a single  multiposition 
switch  in  the  DAASM  system.  Depending  on  the  selected  sub-array,  the  DAASM 
system  samples  each  antenna  sequentially  with  each  transmitted  pulse.  After 
completing  the  sequence  of  G or  12  antennas,  the  cycle  is  repeated  until  512 
samples  are  acquired  on  each  antenna  of  the  selected  array.  In  fact,  the  DAASM 
system  operates  at  any  two  selected  radio  frequencies  at  the  same  time. 

2.  Richard,  D.  W.  (1972)  Twenty- Element  Receive  Array  for  the  DAASM  Ex- 

periment.  Instrumentation  Paper  No.  178. 

3.  Bibl,  K.  (1973)  Doppler /Angle  of  Arrival  Spectral  Measurement  System, 

A FCRL-T  R - 7 3-  075  9 . 
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time  and  space.  In  order  to  accomplish  this,  an  array  of  small  receiving  antennas 
was  erected  in  a linear  configuration  broadside  to  the  required  receiving  direction. 
In  this  case  the  boresight  of  the  array  at  Goose  Ray,  that  is,  the  direction  normal 
to  the  length  of  the  array,  was  magnetic  North,  328°  T.  In  order  to  provide  maxi- 
mum flexibility  in  terms  of  spatial  sampling,  a complex  antenna  arrangement  was 

2 3 

developed  using  a 20-element  array.  The  DAASM  electronic  system  was  de- 
signed to  use  subsets  of  this  array,  consisting  of  either  12  or  6 elements.  Each 

2 

subset  has  either  a uniform  spacing  or  a Minimum  Redundancy  Array  (MRA). 

The  MRA  provides  all  spacings  that  are  multiples  of  a basic  spacing  up  to  some 
maximum  spacing  with  the  minimum  number  of  antennas.  In  the  DAASM  system 
the  several  antenna  sub-array  configurations  available  are  briefly  described  in 
the  following  table: 


Table  1.  DAASM  Arrays  (Dargel 


No.  of 
Elements 

Min.  Spacing 

Max.  Spacing 

Remarks 

Configuration 

Code 

12 

10  m 

110m 

Uniform 

5/0 

1 1 

10  m 

450  m 

M RA 

5/4 

The  configuration  code  above  is  designated  by  the  symbols  T/A,  where  T 
selects  the  frequency  mode  and  sampling  sequence  and  A selects  the  antenna  se- 
quence. As  will  be  described  later,  the  DAASM  svstem  was  designed  to  operate 
with  pulsed  transmissions  at  a 400  Ilz  repetition  frequency.  The  experiments 
described  in  this  report  all  involve  measurements  carried  out  in  conjunction  with 
the  transmitter  on-board  the  KC-135  aircraft  as  a signal  source  flying  to  the 
north  of  Goose  Ray.  Recause  of  transmitter  average  power  limitations  on  the 
aircraft,  pulse  repetition  frequency  was  limited  to  200  Hz,  This  necessitated  a 
change  in  the  Goose  Ray  antenna  system.  Five  additional  sub-arrays,  each  with 
six  elements,  were  devised.  Table  2 describes  these  sub-arrays.  All  of  the 
arrays  are  included  in  the  20  antennas  and  are  selectable  by  a single  multiposition 
switch  in  the  DAASM  system.  Depending  on  the  selected  sub-array,  the  DAASM 
system  samples  each  antenna  sequentially  with  each  transmitted  pulse.  After 
completing  the  sequence  of  6 or  12  antennas,  the  cycle  is  repeated  until  512 
samples  are  acquired  on  each  antenna  of  the  selected  array.  In  fact,  the  DAASM 
system  operates  at  any  two  selected  radio  frequencies  at  the  same  time. 

2.  Richard,  D.  W.  (1972)  Twenty-Element  Receive  Array  for  the  DAASM  Ex- 

periment. Instrumentation  Paper  No.  178. 

3.  Ribl,  K.  (1973)  Doppler/Angle  of  Arrival  Spectral  Measurement  Svstem, 

AFCRL-T  R-7  3-0759. 
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Table  2.  Additional  Sub-Arrays 


No.  of 
Elements 

Min.  Spacing 

Max.  Spacing 

Remarks 

Configuration 

Code 

6 

10  m 

50  m 

Uniform 

7/0 

6 

20  m 

100  m 

Uniform 

7/2 

6 

20  m 

320  m 

M RA  (Approx) 

7/4 

6 

30  m 

720  m 

MRA  (Approx) 

7/6 

6 

70  m 

910  m 

MRA  (Approx) 

7/7 

1 he  actual  procedure  is  to  cycle  through  the  array  the  first  time  using  the  first 
frequency,  then  a second  time  through  the  array  using  the  second  frequency.  On 
the  third  sampling  sequence  through  the  array,  the  first  frequency  is  used  again 
and  this  cycling  is  repeated  until  the  b 12  time  samples  at  each  antenna  and  at  each 
frequency  have  been  recorded. 

Ihis  technique  results  in  a sampling  rate  of  approximately  16  Hz  per  antenna, 
per  frequency,  and  a 32  sec  time  sample  when  the  system  prf  is  200  Hz. 

The  signals  are  recorded  synchronously  on  digital  magnetic  tape  from  four 
sampling  gates  delayed  in  time  from  the  transmitter  pulse  by  a selectable  amount. 
The  output  of  the  four  gates  is  digitized  with  seven  bits  and  then  recorded.  Each 
sample  is  recorded  in  quadrature  at  the  100  kHz  IF.  For  the  experiments  des- 
cribed here,  the  delayed  gates  are  set  at  the  time  of  the  arriving  signals  from  the 
aircraft  which  is  equipped  with  a similar  DAASM  system  and  operates  synchron- 
ously with  Goose  Bay. 

The  spatial  (antenna  array)  and  temporal  samples  form  the  two-dimensional 
array  which  is  used  in  the  data  analysis  to  generate  the  two-dimensional  DAASM 
maps  as  described  in  detail  in  the  next  section. 

The  experiment  described  in  this  report  consists  of  six  aircraft  missions 
covering  the  period  of  December  1973  to  .July  1974.  Each  flight  was  approxi- 
mately 8 hours  in  duration,  and  several  different  flight  patterns  were  executed 
flying  out  to  distances  of  3000  km  to  the  north  of  Goose  Bay.  The  details  of 
these  plans  will  follow  in  Section  2.3  . Whenever  obliquely  propagated  signals 
from  the  aircraft  were  available,  the  delayed  sampling  gates  were  made  to  coin- 
cide with  the  received  signal  and  the  data  digitized  and  recorded.  Two  frequen- 
cies were  selected  for  each  half  hour  of  the  flight  before  take-off  using  the  best 
propagation  predictions  available.  Jr  general,  each  predicted  frequency  was 
used  for  at  least  1 hour.  With  the  availability  of  real-time  oblique  ionograms 
(Section  2.2)  between  aircraft  and  Goose  Bay,  the  selection  of  operating  frequency 
was  often  made  during  the  flight. 
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2.2  Oblique  lonograms 

In  November  1973,  it  was  decided  to  generate  oblique  stepped  frequency  iono  • 
grams  between  the  KC-135  aircraft  and  the  Goose  Ray  Ionospheric  Observatory 
during  the  same  time  that  the  routine  backscatter  ionogram  was  being  made  at 
Goose  Bay.  Roth  the  Goose  Bay  and  the  aircraft  sounding  systems  begin  the 
stepped  frequency  mode  simultaneously  and  each  records  an  oblique  ionogram. 

The  frequency  range  is  normally  6 to  16  MHz  in  100  kHz  steps  and  the  ionogram 
takes  approximately  1 min  (Figure  1).  These  ionograms,  although  instituted  pri- 
marily for  the  purpose  of  studying  the  propagation  in  the  polar  regions,  have 
played  a major  role  in  the  DA  ASM  experiment.  First  they  have  made  possible 
the  real-time  selection  of  the  two  operating  frequencies.  By  agreement  with  the 
operating  personnel  of  both  the  ground  station  and  the  aircraft,  and  after  inspection 


‘ ■"  ' 1 5 s-  5 1.  ..5.  ..  .5  ..s 

I 1 1 1 1 1 I I I I I 

6'l  7.1  8-1  9.1  10.1  III  12.1  13.1  14.1  15.1  16.1 

FREQUENCY  IN  MHz 

Figure  1.  Oblique  Ionogram,  Goose  Bay, 

Flight  4-196,  15  .1  ul  1974:  1641  ITT:  1240  km 
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of  the  oblique  ionogram,  certain  specified  changes  are  allowed  at  one  end  of  the 
path  and  the  operator  at  the  other  end  attempts  to  follow  these  changes.  With  the 
two-way  oblique  ionograms.  this  real-time  method  has  worked  exceedingly  well 
and  has  resulted  in  frequency  usage  closer  to  the  optimum  than  is  possible  using 
4 

predictive  methods. 

The  second  benefit  derived  from  the  oblique  ionogram  has  been  in  the  post 
analysis,  specifically  the  possibility  of  making  mode  identifications.  An  important 
part  of  the  DAASM  program  is  to  relate  the  observed  spectral  character  of  the 
fixed  frequency  maps  to  the  propagation  mode.  It  is  through  understanding  this 
relationship  that  techniques  for  improving  OTH  radar  performance  can  be  de- 
veloped. 5 

The  oblique  ionograms  are  made  using  the  existing  Digisonde  128  systems 

located  at  Goose  Ray  and  in  the  AFCRL  KC-135  aircraft.  The  data  is  recorded 
on  digital  magnetic  tape  as  well  as  displayed  on  site  in  real  time  in  the  manner 
used  for  all  "Digisonde"  data. 


2.3  Flight  Flans 

For  each  of  the  six  flights  used  for  this  report,  a flight  plan  was  developed 
for  use  by  the  aircraft  navigator.  The  actual  flight  paths  deviated  only  slightly 
from  the  planned  path,  and  in  the  subsequent  description  only  the  actual  paths  are 
presented.  The  chosen  paths  were  often  a compromise  between  competing  objec- 
tives Certain  compromises  in  Right  path  were  sometimes  made  to  permit  auro- 
ral oval  monitoring  by  other  scientists  involved  in  the  airborne  measurements 

program. 

Each  of  the  six  Right  paths  are  presented  in  three  different  manners.  ie 
first  display  for  each  flight  shows  the  actual  geographic  location  of  the  aircraft  as 
a function  of  time  on  a geomagnetic  projection,  while  the  second  display  shows 
range  and  bearing  of  the  aircraft  from  Goose  Bay.  The  third  display  shows  the 
path  of  the  aircraft  in  corrected  geomagnetic  time  vs  latitude.  The  auroral  oval 
for  Q=  3 remains  fixed  in  this  coordinate  system.  This  allows  the  user  to  see 
the  physical  relationship  between  the  aircraft.  Goose  Bay  and  the  oval  during  the 
Right.  These  aircraft  track  maps  are  useful  when  interpreting  the  observed 
propagation  modes  and  DAASM  maps. 


4 Barghausen,  A.  F. . et  al  (196'91  Predicting  Long  Term  Operational  Pa  ram  - 
' eters  of  High  Frenuencv  Sky-Wave  Telecommunications  Systems.  RsSA 
Tech.  Report  ERL  llO-ITS-'/B. 

5.  Bibl,  K..  et  al  (1970)  Digital  Interpreting  Goniometric  Ionospheric  Sounde£, 
AFCRL- 7 1-0002. 
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The  six  flights  used  in  this  report  are  detailed  in  the  following  table: 

Table  3.  DA  ASM  Flights 


DAASM  No. 

Date 

Take-Off  (UT) 

Landing  (UT) 

4- 

1 

O 

CJl 

27  Feb  1974 

1900 

034  2 

DM  4-060 

1 Alar  1974 

2355 

0829 

DAI  4-141 

21  Alay  1974 

1743 

0220 

DM  4-142 

22  Alay  1974 

2032 

0457 

DM  4-193 

1 2 .1  ul  1974 

1440 

2250 

DM  4-196 

15  .1  ul  1974 

1440 

2236 

3.  ANALYTIC  TKClIMOl T.S 
3.1  DA  ASM 

■As  stated  m the  introduction,  the  purpose  of  the  DA  ASM  project  is  to  provide 

# deSCripti°n  °f  ,he  Character  °f  r;idio  «‘*na1s  arriving  from  a point  source  through 
a relatively  disturbed  ionospheric  medium.  In  order  to  provide  such  a description 
the  approach  chosen  is  to  devise  a signal  power  density  estimate  as  a function  of  two 
enables.  Doppler  frequency,  <o.  and  wave  number,  k.  associated  with  arrival  angle. 

1 he  power  as  a function  of  these  two  variables  results  from  the  two- 
dimensional  transforms  of  the  recorded  amplitude  data  in  space  and  time,  that  is, 
x and  t.  The  DAASM  system  at  Goose  Hay  utilizes  a one-dimensional  linear 
azimuthal  arrival  angle.  The  data  is  ambiguous  with  respect  to  the  elevation 
angle  of  the  source  or  sources. 

I hree  methods  of  estimating  the  power  density  function  are  described.  T le 
three  methods  arose  from  the  various  antenna  configurations  available  with  the 
DAASM  system.  The  third  method  is  the  only  one  used  for  all  th'  data  presented 
here,  even  though  it  is  relatively  time  consuming  on  the  computer  because  it  in- 
volves matrix  inversion. 

As  a general  introduction  to  the  analytic  techniques  used  here,  a short  des- 
cription of  the  basic  approach  is  given. 

If  the  time  samples  are  understood  to  constitute  one  dimension  and  the  spa- 
tial samples  to  be  the  second  dimension,  then  the  two-dimensional  Fourier  trans- 
forms yields  the  Doppler  frequency/wavenumber  power  density  (FWPD)  spectrum. 

Thus,  if  A.(tu)  and  AJu)  denote  the  discrete  complex  Fourier  spectral  line 
at  Doppler  frequency,  co  at  antennas  j and  m respectively,  then  the  cross  spec- 
trum between  the  antennas  is  given  by: 
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S (CO)  = A.(co)  • A'  (co)  . 
j m j m 


In  general,  a smoothed  version  of  cross-spectral  estimates  S.  (co)  is  always 
used.  The  estimated  FWPD  spectrum  is  defined  by 


- V 

P(co.k)  = Z 
j=l 


N 

z 

m=  1 


W.W  S.  (CO) 
J m jm 


exp  i k • (r.-r  ) 

J m 


(1) 


where  N is  the  total  number  of  antennas,  ?n's  are  the  vector  distances  from 


some  arbitrary  origin,  k is  the  wavenumber  vector,  and  the  W's  are  suitable 

.th 


weights  imposed  on  the  il“  antenna. 
1 


In  a report,  different  methods  of  spectral  analysis  to  extract  Doppler  fre- 
quency, phase  and  coherence  information  are  described  in  detail.  It  was  found 
that  the  frequency-averaging  technique  provided  a very  efficient  and  reliable 
means  of  cross-spectral  estimation,  and  only  the  frequency-averaging  technique 
of  estimating  the  cross-spectra  is  used  in  all  of  the  data  involved  here.  Once  the 
cross-spectra  have  been  estimated,  the  spatial  component  of  the  FWPD  spectrum 
can  be  estimated  by  several  methods.  The  three  methods  of  wavenumber  analysis 
can  be  delineated  in  terms  of  the  choice  of  the  weights  W.  and  are:  (1)  The  Fast 
Fourier  Frequency-Wavenumber  Analysis  or  the  Conventional  Analysis,  which  is 
most  useful  for  an  initial  diagnostic  scan  in  the  frequency-wavenumber  (co- k)  space. 
Wavenumber  sidelobe  characteristics  can  be  varied  by  proper  choice  of  weights 
for  the  different  cross-spectral  terms;  however,  the  resolution  is  relatively  poor. 

(2)  The  minimum  redundancy  version  of  the  conventional  method  is  only  applicable 
to  the  data  from  a strict  minimum  redundancy  configuration  of  the  antennas, 

(3)  The  adaptive  method  using  the  maximum  likelihood  principle  offers  high  reso- 
lution in  spatial  characteristics  and  is  applicable  to  any  arbitrary  antenna  config- 
uration. Unlike  the  other  two  methods,  this  method  assigns  adaptive  weights  to 
the  different  cross-spectral  terms  depending  on  the  azimuthal  look-angle  and  the 
spatial  characteristics  of  interference  and  noise. 

The  part  of  the  analysis  and  computation  common  to  all  three  methods  is  the 
determination  of  a suitable  cross-spectral  estimate  between  any  given  pair  of 
antennas. 

The  recorded  data,  512  points  per  subcase  for  each  antenna,  consists  of  two 


sets  of  256  equally-spaced  complex-valued  time  samples  displaced  from  each 
other  by  a ratio  (3  = • bach  256  sample  set  could  be  used  for  processing;  how- 

ever, it  is  practical  to  treat  only  the  middle  128  samples  in  each  set  in  order  to 


Pfister,  W.  , et  al  (1975)  Pulse  Sounding  with  Closely  Spaced  Receivers  as  a 
Tool  for  Measuring  Atm^i^7inc  Motions  and  Fine  Structure  in  the  Iono- 
sphere, Vol.  VII,  Environmental  Research  Paper  Nlo.  506,  April  1975. 
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allow  for  adjustment  of  range  gates  and  frequency  before  processing  the  data. 

Thus,  the  sets  consisting  of  128  samples  in  each  subcase  are  immediately  trans- 
formed into  the  frequency  domain  via  a Fast  Fourier  Transform.  The  resulting 
128  complex-valued  discrete  Fourier  spectrum  spans  a basic  Doppler  frequency 
range  of  8 Hz.  It  is  possible  to  either  combine  or  interlace  the  resulting  two  sets 
of  Fourier  spectra.  If  the  Fourier  spectra  are  aligned,  the  signal-to-noise  ratio 
can  be  improved  but  the  basic  Doppler  frequency  range  is  still  8 Hz.  If  the  two 
sets  are  interlaced,  the  Doppler  frequency  range  can  be  doubled  to  1G  Hz  without 
any  improvement  in  the  signal-to-noise  ratio. 

In  tne  case  of  signals  transmitted  bv  the  aircraft,  the  basic  range  of  8 Hz  is 
sufficient.  Hackscatter  signals  often  occupy  a band  wider  than  8 Hz,  in  which  case 
interlacing  proves  to  be  very  useful.  In  any  case,  the  combined  or  interlaced 
Fourier  spectra  form  the  starting  point  for  further  data  processing. 

Three  types  of  output  have  been  designed.  These  consist  of  (1)  power  spectra, 
(21  coherence  spectra,  and  (31  maps  of  Doppler  frequency  vs  arrival  angle  called 
DAASM  maps. 

3.1.1  POWl-IR  SPKCTRA 

Power  spectra  are  calculated  for  each  antenna  by  taking  the  product  of  each 
Fourier  spectral  line  with  its  own  complex  conjugate.  Of  course,  the  Fourier 
spectra  are  normalized  in  each  channel  so  that  the  total  power  in  each  channel 
equals  unity.  That  is,  if  g^  represents  the  unnormalized  mth  frequency  line, 
then  the  normalized  Fourier  spectrum  is  simply 


The  power  spectra  are  smoothed  by  the  frequency-averaging  technique'  over 
seven  frequencies  concurrently  with  the  operation  of  evaluating  the  product  of  the 
Fourier  spectra.  This  feature  enables  the  cross-spectra  to  be  estimated  in  an 
identical  way  merely  by  inserting  the  Fourier  spectra  of  the  appropriate  pair  of 
antennas.  The  ffower  spectra  and  the  cross-spectra  are  different  only  in  a mat- 
ter of  the  antennas  involved. 

The  power  spectra  thus  obtained  are  pic  tted  on  a standard  format  of  logarith- 
mically-scaled numbers  ranging  from  0 to  15  which  corresponds  to  a linear 
scale  from  0 to  HP.  The  exact  equation  defining  the  scales  is 

1.  Pfister,  W.  , et  al  (19751  Pulse  Sounding  with  Closely  Spaced  Receivers  as  a 
Tool  for  Measuring  Atmospheric  Motions  and  Fine  Structure  in  the  lono- 
sphere,  Vol.  V(l,  Environmental  Research  Paper  No.  506,  April  1975. 
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I he  power  spectral  outputs  provide  the  basis  for  selection  of  data  segments 
suitable  for  the  more  complex  and  time  consuming  P(co.  k)  analysis.  The  data 
consists  of  subcases  where  each  subcase  has  four  different  range  gates  and  two 
carrier  frequencies.  Therefore,  for  the  six  antenna  configurations  there  will  be 
4!i  channels  of  power  spectra  for  each  subcase. 

:i.i.2  ro!ii;m:\ri: 

Coherence  spectra  between  the  mth  and  the  nth  channels  (antenna)  is  defined 
by  the  equation: 


fronted)  I = 

m n j 


lSrnnU0)i 
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WllCre  Smn  ,aJ)  is  ,h0  -smoothed  «-ross  spectrum  at  Dormler  frequency  w between  mth 

:‘nd  " Ch:mnels  ;‘nd  wo“ld  '*  the  power  spectrum  of  the  mth  cl  annel.  It 

is  clear  that  the  coherence  is  ..  number  which  lies  between  0 and  1 . This  linear 
range  is  converted  to  an  inverse  hyperbolic  tangent  scale  and  multiplied  bv  a suit- 
able scale  factor  to  provide  the  printed  output  numbers  for  the  data  in  Section  4 
The  following  table  shows  the  relationship  between  coherence  value  and  output 
numbers.  The  coherence  spectra  are  plotted  for  each  available  distance  between 
the  various  pa,  s of  antennas  m order  of  increasing  distance.  Since  the  coherence 
is  only  available  at  specified  distances,  the  plots  are  backfilled  at  intermediate 
distances  with  the  same  value.  It  was  found  that  the  raw  coherence  thus  plotted 
appeared  sometimes  to  be  a discontinuous  function  of  distance  instead  of  a mono- 
ton,cally  decreasing  function.  Therefore,  three  successive  subcases  have  been 
averaged  lor  each  distance  at  each  Doppler  frequency  to  provide  a much  smoother 
Picture  and  thus  improve  the  statistical  confidence  limits  which  depend  on  the 

sample  size.  A least-square  linear  fit  is  plotted  for  a sample  of  these  data  in 
I*  igure  2. 


The  coherence  spectra  are  primarily  useful  in  determining  the  maximum  use 
ful  aperture  as  a function  of  operating  frequency  and  the  signal  propagation  mode. 
3.  1.  3 DA  ASM  MAPS 

DA  ASM  maps  consist  of  plots  of  power  density  vs  Doppler  frequency  and 
wavenumber  which,  in  turn,  is  simply  related  to  the  azimuthal  angle  of  arrival. 

3.  1.3.  1 Bartlett  Estimate 


A particularly  efficient  algorithm  to  calculate  the  (go- k)  spectra  is  achieved 
if  the  weights  in  Eq.  (1)  are  chosen  to  correspond  to  a triangular  taper  function 
anc  will  be  referred  to  as  Fast  Frequency- Wavenumber  analysis  or  the  general- 
ized Bartlett  estimate. 

In  such  a case,  the  individual  cross-spectral  terms  need  not  be  separately 
calculated.  The  separate  Fourier  spectra  from  each  antenna  output  can  be  used 
directly,  and  frequency  smoothing  can  be  applied  in  the  final  step.  The  details  of 
the  algorithm  are  as  follows: 


i Aj(co)  exp  i 0.(w)  represent  the  Fourier  spectral  line  at  frequency  go 


from  the  jth  antenna,  then  the  cross-spectrum  between  the  ith 


and  m antennas 


S. 

im 


(go) 


A (go)  | I A (go)  I exp  i (0.  - 0 ) 
I m i m 


and  the  (go-  k)  spectrum 
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P(co,k)  = X'  S.  (go)  exp 
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where  k is  the  wavenumber. 
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Figure  2.  Coherence  vs  Distance,  lypical  Plot 
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Eq.(2)  can  be  written  alternately  as 
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(3) 


The  subscripts  j and  m can  be  replaced  by  a single  subscript  n m-j , and 
the  double  summation  can  be  replaced  by  a singlfe  summation.  Thus, 


P(co,k)  = 


A (oj  > exp  i 0 (oj  ) - k 
n n 


(41 


Eq.  (41  can  be  computed  much  faster  than  Eq.  (31. 

— > — * 

P(o),  kt  is  calculated  at  as  many  points  as  desired  in  k,  say  with  increments 

•*> 

Ak.  Each  successive  point  in  k can  be  computed  from  the  previously  computed 
point  by  simple  complex  multiplication  because  exp  i (k  4 A.':)  = exp  ik  • exp  i Ak  . 
This  step  is  common  to  all  three  approaches  of  (io-k)  spectral  calculations  and 
saves  considerable  computer  memory  and  time. 

Since  frequency  smoothing  is  independent  of  the  operations  of  k , a set  of 
weights  \V(u>)  is  calculated  using  P((U,o)  for  each  co,  and  each  P(aj,k)  is  multi- 
plied by  the  appropriate  W(a>)  to  yield  P(w,k),  the  frequency  smoothed  estimate. 
The  relation  between  the  wavenumber  k and  the  direction  of  arrival  0 measured 
from  boresight  for  the  case  of  a linear  array  is  given  by 


k = 


k 


2 TTf 

c 


sin  0 


where  k is  now  a scalar,  f the  operating  carrier  frequency,  and  c the  velocity 
of  light. 

3.  1.3.2  Hanning  Estimate  (Minimum  Redundancy  Version) 

The  minimum  redundancy  array  version  can  be  treated  as  a special  case  of 
the  general  equation  for  the  (w-k)  spectra  given  by  Eq.  (1).  It  is. assumed,  of 
course,  that  a certain  number  of  consecutive  multiples  of  the  basic  antenna  spacing 
are  available  by  suitably  placing  the  different  antennas.  Distances  beyond  the 
largest  consecutive  multiple  are  considered  as  having  zero  weight;  that  is,  their 
contribution  to  the  (co-k)  spectra  is  zero.  The  non-zero  weights  decrease  as  a co- 
sine squared  function  with  distance.  In  addition,  the  self-spectral  term  is  assigned 
a weight  of  0.5,  so  that  the  expression  for  the  (w-k)  spectrum  can  be  written  as 
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•‘'here  Sq  corresponds  to  the  self  term  anti  S;  is  the  cross -spectra]  term  corre- 
sponding to  the  pair  of  antennas  separated  by  a distance  d^  . L is  the  largest  suc- 
cessive multiple  of  the  basic  spacing  and  W.  is  the  weight  given  by  the  cosine  - 
squared  function.  It  is  well  known  that  ' he  cosine -squared  weights  result  in  good 
sidelobe  suppression.  It  can  be  shown"*  that  for  6 antennas  only  13  consecutive 
multiples  are  possible,  and  for  1 1 antennas  the  corresponding  number  is  45.  A 
similar  configuration  for  more  than  11  antennas  is  not  known  at  the  present  time, 
hach  spacing  occurs  only  once,  and  hence  the  name  "minimum  redundancy"  for  the 
array.  Only  those  cross -spectral  terms  corresponding  to  the  multiples  of  the 
basic  spacing  need  be  calculated.  Thus,  for  1 1 antennas,  45  of  66  possible  terms 
are  available.  The  applied  weights  are  almost  optimal  for  the  configuration  and 
there  is  excellent  sidelobe  suppression,  if  the  carrier  frequency  is  chosen  such 
that  the  basic  spacing  is  a half-wave  length  or  smaller,  then  there  will  be  no  grat- 
ing lobes. 

3.  1.3.3  Maximum  Likelihood  Kstimate 

I he  third  approach  to  (cu-k)  spectral  estimation  employs  an  adaptive  scheme 
to  generate  the  set  of  weights  in  the  filter  given  by  Eq.  (1)  Instead  of  using  a 
fixed  set  of  weights  for  all  k , the  adaptive  scheme  allows  fcr  the  set  of  weights  to 
be  varied  depending  on  the  particular  k as  well  as  the  character  of  the  signal, 
that  is,  interference  and  noise.  In  other  words,  the  set  of  weights  is  dictated  by 
the  data  sample  available,  subject  to  the  constraints  that  a unit  plane  wave  signal 
travelling  with  a given  wavenumber  k = kQ  (arriving  from  the  direction  towards 
which  the  antenna  beam  is  pointed),  is  passed  undisturbed  through  the  spatial  filter, 
but  everything  else  for  which  k = kQ  (arriving  from  other  directions)  is  rejected  in 
an  optimal  manner.  Thus,  the  true  spectrum  is  reproduced  with  high  resolution 
and  not  smeared,  which  normally  occurs  as  a result  of  the  convolution  between 
the  true  spectrum  and  the  particular  array  pattern.  This  choice  of  constraints  re- 
sults in  a convenient  mathematical  formulation.  It  is  only  necessary  to  calculate 
the  inverse  of  the  matrix  of  cross-spectral  terms  for  each  Doppler  frequency  of 
interest  in  order  to  furnish  the  optimal  set  of  weights  automatically  adapted  to 
each  k for  a given  data  sample.  The  computational  procedure  is  set  forth  in  the 
following  steps. 

The  normalised  input  Fourier  spectra  from  the  various  antennas  is  converted 
to  cross -spectral  estimates  .8.^(0))  by  any  one  of  several  spectral  estimation  pro- 
cedures, in  this  case  by  frequency  averaging  and  smoothing  with  seven  frequencies. 

6.  Moffet,  A.T.  (1968)  IEEE  Transactions  on  Antennas  and  Propagation  AP-  16 
(No.  2>:  172.  — “ 


Subscripts  j and  m refer  to  the  antennas  in  question.  Next,  an  NXN  Matrix  of 
cross -spectral  estimates  fS(o!)]  is  organized  at  each  to,  where  N is  the  total  num- 


ber of  antennas.  Not  all  NXN  terms  need  be  evaluated  since  S.  (oo)  = S*  .(to). 

jm  mj  ’ 

wdiere  the  asterisk  denotes  complex  conjugation.  The  sum  of  the  diagonal  elements 


in  the  matrix  is  stored  as  a measure  of  total  power  at  anv  given  Doppler  frequency 

and  used  later  for  renormalization.  The  matrix  itself  is  normalized  by  replacing 

1/2 

each  element  S.^w)  X Sjm  co)/[S..{W)  S^Co)]  , a step  which  ensures  sensor 
equalization.  The  normalized  cross -spectral  matrix  for  Doppler  frequency  to  is 
then  inverted.  The  elements  of  the  inverted  matrix  are  used  in  the  equation  for 
(to-k)  spectra 
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It  is  seen  that  the  weights  VV.  and  the  cross -spectral  elements  [Eq.(l)]  together 
ace  replaced  by  the  corresponding  inverted  matrix  elements  and  the  reciprocal  is 
taken  to  yield  the  high-resolution  (co-k)  power  density  estimate.  The  increments 
over  k are  calculated  in  a chain  sequence  from  the  various  values  in  the  same 
manner  as  indicated  before.  The  final  step  is  to  restore  relative  power  levels  at 
the  different  Doppler  frequencies.  First,  the  values  of  P(co,k)  are  summed  over 
all  k and  the  normalized  values  of  P(ce,k)  are  obtained  by  simply  dividing  by  this 
sum.  Next,  the  normalized  P(u,  k)  are  multiplied  by  the  sum  of  the  diagonal  ele- 
ments of  the  original  cross -spectral  matrix  which  have  been  stored  as  indicated 
before.  I hus,  the  renormalized  P(to,k)  are  obtained  with  an  appropriate  scale 
factor  and  are  recorded  on  magnetic  tape  for  later  printouts  (see  Section  4.2.3). 


3.2  Oblique  lonoerams 

During  the  six  flights,  the  oblique  ionograms  have  been  scaled  for  both  mode 
availability  and  frequency.  The  oblique  ionograms  are  taken,  ore  every  5 min, 
which  represents  a considerable  effort.  The  results  are  used  to  provide  the  iden- 
tification necessary  to  relate  the  DAASIV1  maps  with  propagation  mode  Another 
aid  to  mode  identification  is  the  Polar  Model  ITS-78  program.7  The  polar  model 
program  has  been  modified  for  use  with  aircraft  to  Goose  Bay  propagation  paths, 
that  is,  where  one  end  of  the  propagation  path  is  continuously  changing.  A syn- 
thesized ionogram  is  calculated  from  this  program  and  compared  with  the  measured 
one. 


The  calculated  and  measured  oblique  ionograms  together  provide  a reasonably 
reliable  method  of  mode  identification. 


7.  t,lkins,  r.  l.  (1973'  An  Empirical  Model  of  the  Polar  Ionosphere.  Survev  in 
Geophysics,  No.  267.  c 
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4.  RESULTS 


4.1  General 


I he  results  of  the  aircraft  measurements  will  be  presented  in  sequential  sec 
tions  with  each  flight  as  a complete  package 

data  will  be  presented  in  a subsequent  report,  and  here  only  data  w 
to  illustrate  the  character  of  the  aircraft  signals  propagated 
Within  each  flight  section 
ures  3(a) 


The  detailed  analysis  of  the  DA  ASM 
ill  be  presented 
over  auroral  paths, 
the  different  types  of  display  will  be  indicated.  Fig- 
lb),  (c)  through  8(a),  (b),  (c)  for  Flights  4-058,  4-060,  4-141,  4-142, 
4-193,  ami  4-196  are  concerned  with  the  aircraft  track  and  propagation  paths. 

A limited  quantity  of  data  is  presented  here,  varying  from  eight  maps  for 
Flight  4-196  to  two  maps  for  Flight  4-14  1.  These  maps  are  given  in  Figures  9(a) 
and  (b)  through  14(a)  and  (b).  The  maps  were  selected  to  illustrate  certain  nenera 


(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  3.  DA  ASM  Flight.  4-058,  27  Feb  1974 
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AIRCRAFT  RANGE  (Km) 


(b)  Range  and  Bearing  vs  Time  of  Aircraft  from  Goose  Bay 


1 2 

'0 

70' 


(c)  Aircraft  Tracks  Relative  to  Auroral  Oval  (Q  = 3),  in 
Geomagnetic  Bocal  Time  and  Latitude 

Figure  3.  DA  ASM  Flight  4-05B,  27  Feb  197  4 
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AIRCRAFT  BEARING 


features  of  the  observed  behavior  and  special  events.  For  each  time  selected,  the 
maps  are  included  for  both  Doppler/Arrival  Angle  and  Doppler  Frequency/Coher- 
ence Distance.  This  results  in  a rather  complete  description  of  the  received  sig- 
nals, The  details  of  the  format  for  this  map  data  will  be  given  in  Sections  4.2,3 
and  4.2.4.  For  each  flight,  the  maps  will  begin  with  Figure  9(a)  for  Doppler/ 
Arrival  Angle  and  Figure  9(b)  for  Doppler/Coherence  Distance  maps.  Each  figure 
will  consist  of  all  the  maps  for  that  flight  identified  by  time  and  arranged  sequentially. 

The  map  data  for  each  I light  is  supplemented  by  the  geographical  and  auroral  oval 
plots  of  the  aircraft  track.  The  auroral  oval  aircraft  track  shows  the  relative  posi- 
tion of  the  aircraft,  Goose  Bay,  and  the  auroral  oval  at  all  times  during  the  flight. 

4.2  Data  Format 

4.2.  1 AIRCRAFT  TRACK  - OEOGRAPFirCA  I.  COORDINATES 

Using  the  Kf  -135  aircraft  navigator  log,  the  position  of  the  aircraft  at  each 
observation  time  is  plotted  in  geographic  coordinates  in  F'igures  3(a)  through  8(a) 
for  each  flight  on  a polar  projection  of  the  Earth.  Along  the  track  are  indicated  the 
hourly  positions.  In  order  to  make  some  of  the  more  complicated  flight  paths  clear 
and  for  use  in  analysis  of  the  DAASM  maps.  Figure  2(b)  is  n plot  of  range  and 
bearing  of  the  aircraft  measured  from  Goose  Bay  as  a function  of  time. 


(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  4.  DA  ASM  Flight  4-060,  1 Mar  1974 
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(b)  Range  and  Bearing  vs  Time  of  Aircraft  from  Goose  Bay 
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(c)  Aircraft  Tracks  Relative  to  Auroral  Oval  (Q  = 3),  in 
Geomagnetic  Local  Time  and  Latitude 
Figure  4.  I)A ASM  Flight  4-OCO,  1 Mar  1974 
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4.  2.  2 AIRCRAFT  TRACK  - AURORAL  OVAL  COORDINATES 


In  order  to  see  the  relationship  between  the  aircraft  position  and  the  auroral 
oval  as  a function  time  during  the  flight,  a second  plot  of  aircraft  position  in  geo- 
magnetic latitude  against  corrected  geomagnetic  time,  in  polar  form,  is  shown  in 
Figures  3(c)  through  8(c).  In  this  coordinate  system,  the  auroral  oval  remains 
fixed  while  the  receiving  site  at  Goose  Bay  circles  at  constant  geomagnetic  lati- 
tude, 53.3°N,  and  is  shown  at  hourly  positions.  In  Figure  3(c),  the  one-hop  and 
two-hop  reflection  points  are  indicated  for  the  aircraft  to  Goose  Bay  at  0000  UT. 
The  points  labelled  1 and  3 are  the  two-hop  reflection  points,  and  the  point 
labelled  2 is  the  one-hop  reflection  point.  It  is  obvious  from  Figure  3(c)  that  the 
one-hop  reflection  point  lies  within  the  oval  near  the  poleward  edge  while  for  the 
two-hop  mode,  one  reflection  point  lies  in  the  polar  cap  and  the  second  reflection 
lies  within  the  oval.  This  information  will  be  useful  for  one  interpretation  of  the 
DAASM  maps. 

l . ✓ ^ 11  ^ V, 
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(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  5.  DAASM  Flight  4-141,  21  May  1974 
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(b'  Range  and  Hearing  vs  Time  of  Aircraft  from  Goose  Hay 
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(c'  \ i r craft  Tracks  Relative  to  Auroral  Oval  (Q=  31,  in 

Geomagnetic  Local  Time  and  Latitude 
figure  5.  DA  ASM  flight  4-141,  21  May  1974 
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4.  2.  3 DAASM  Maps 

At  any  selected  time  and  range,  a DAASM  map  is  generated  according  to  the 
technique  described  in  Section  3.  1.  3.  3 . These  maps,  for  each  of  the  two  selected 
operating  radio  frequencies,  represent  the  Doppler  frequency  vs  arrival-angle 
spectrum  of  the  aircraft  signal.  The  two  operating  frequencies  are  displayed  one 
above  the  other,  and  the  abscissa  is  the  Doppler  frequency  in  1 Hz  steps  about  a 
chosen  center  frequency  indicated  in  the  legend  attached  to  each  map.  The  total 
Doppler  frequency  range  is  always  8 Hz  for  the  maps  included  in  this  report.  The 
ordinate  scale  represents  the  arrival  angle,  linear  in  wave  number,  measured 
from  boresight  (0 0 1 to  ± 90°  at  the  end  of  each  map.  The  arrival-angle  scale  is 
repeated  twice;  once  for  each  radio  frequency.  The  first  frequency  listed  corre- 
sponds to  the  lower  map. 

a 

Using  a special  printing  system,  the  presented  maps  have  a dynamic  range 
of  approximately  22  dB  with  the  difference  between  adjacent  numbers  equal  to 


(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  6.  DAASM  Flight  4-142,  22  May  1974 
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(b)  Range  and  Hearing  vs  Time  of  Aircraft  from  (loose  Bay 
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(cl  Aircraft  Tracks  Relative  to  Auroral  Oval  (Q=3) 
Geomagnetic  Local  Time  and  Latitude 

Figure  6.  DAASM  Flight  4-  14  2,  22  May  1974 
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1.5  dB.  The  numbering  system  is  designed  so  that  the  larger  numbers  (numbers 
range  from  0 to  15)  appear  more  dense,  giving  a visual  appearance  of  increasing 
blackness  near  the  peaks  of  the  signal.  For  example.  Figure  9(a)  (UT  2343)  shows 
the  DAASM  map  for  aircraft  signals  at  8.2  and  6.6  MHz  over  a distance  of  1960  km. 
The  DAASM  map  in  this  case  was  generated  using  a rather  small  array  with  an 
aperture  of  50  m (see  Section  2.1).  For  frequency  1,  8.2  MHz,  the  received  power 
is  very  close  to  boresight,  0°,  and  at  a Doppler  frequency  of  -5.4  Hz.  The  raw 
Fourier  spectra  were  smoothed  over  seven  frequencies,  and  the  small  antenna 
aperture  results  in  a minimum  width  of  the  displayed  signal  of  the  order  shown. 

For  other  maps  where  the  large  antenna  arrays  were  used,  narrower  angular 
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(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  7.  DAASM  Flight  4-193,  12  ,lul  1974 


1.  Pfister,  W.  , et  al  (1975)  Pulse  Sounding  with  Closely  Spaced  Receivers  as  a 
Tool  for  Measuring  Atmospheric  Motions  and  Fine  Structure  in  the  lono- 
sphere,  VoTT  VII,  Enivornmental  Research  Paper  No.  506,  April  1975. 


spectra  are  often  observed.  Caution  must  be  exercised  in  interpreting  the  ob- 
served angular  spreading  until  the  specific  antenna  array  is  taken  into  account. 
The  bulk  of  the  data  in  this  report  was  taken  using  the  7/6  configuration  with  an 
aperture  of  720  m allowing  relatively  high  angular  resolution. 

Examples  of  narrow  angular  spectra  using  the  7/6  array  are  shown  in  Fig- 
ure 12(a)  (FT  1040).  Here  the  angular  width  is  as  small  as  ; nv  observed  for  sig- 
nals propagated  in  the  Arctic  ionosphere.  The  Doppler  frequency  spreading  is 
also  minimal.  ( .msistent  with  the  Doppler  frequency  smoothing.  The  relatively 
high  amplitude  at  other  arrival  angles  for  the  same  Doppler  frequency  as  the  air- 
craft signal  results  from  the  sidelobes  generated  by  the  "Maximum  Likelihood  1 
processing  (Section  2, 1,2,3),  In  this  case,  the  sidelobes  are  some  la  dll  ht.ow 
the  peak.  This  is  acceptable  considering  the  relatively  few  antennas  (six)  used  in 
these  samples.  The  sidelohe  level  suppression  depends  critically  on  the  number 


(a)  Aircraft  Path  in  Geographic  Coordinates  on  a Polar  Geomagnetic  Projection 
Figure  8.  DA  ASM  Flight  4 - 1 PC . 16  .1  ul  1 074 
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of  antennas  used  in  the  analysis.  All  data  processed  for  this  report  used  six 
antenna  arrays.  The  resulting  high  sidelobe  level  often  makes  it  difficult  to  iden- 
tify the  aircraft  signal;  however,  the  MI.M  always  results  in  t^e  largest  power 
(highest  number'  at  the  correct  signal. 


4.  2.4  CnilKHKWK  MAPS 

The  coherence  between  pairs  of  antennas  within  the  array  is  defined  as  the 


magnitude  of  the  cross-spectra  (Section  .‘1.1.21.  For  each  Doppler  frequency  line 
in  the  spectra  the  coherence  is  calculated  for  all  the  pairs  of  antennas,  that  is, 
all  spacinus  available  in  the  selected  antenna  array.  The  coherence  map  (Fig- 
ures (Kb)  through  14(b!  displays  the  coherence  level  plotted  for  each  Doppler  fre- 
quency vs  distance.  Kach  distance  line  in  the  coherence  map  represents  10  meters 
spacing  up  to  a maximum  of  060  meters.  The  maximum  distance  and  the  available 
intermediate  distances  depends  again  on  the  antenna  configuration  chosen. 

For  example,  the  commonly  used  7/6  antenna  array  uses  six  antennas  with 
physical  spacings  of  20,  60,  100,  200,  and  220  meters  giving  16  possible  spacings 
of  20,  60,  100,  160,  200,  220,  260,  290,  220,  260,  200,  420,  400,  600,  and  720 
meters.  All  of  these  spacings  are  used  for  coherency  calculations  and  for  plotting 
purposes.  The  coherence  for  distances  less  than  20  meters  is  set  equal  to  C Oil 
(20  m)  and  the  coherence  for  distances  of  40  and  50  meters  is  set  eaual  to  COM 
(60  ml.  This  back-filling  continues  up  to  the  largest  available  distance;  in  this 
case  720  meters,  lletween  720  meters  and  060  meters,  the  end  of  the  graph,  no 
data  are  available.  For  the  filled  arrays,  configuration  5/0,  7/0,  and  7/2,  the 
largest  spacing  is  in  general  less  than  100  meters  and  the  coherence  data  of  little 
value,  though  these  coherence  maps  are  included  for  completeness. 

Two  typical  examples  are  discussed  here  for  general  information.  The  first 
example  chosen  illustrates  high  coherence  across  the  entire  large  array  and  the 
second  coherence  which  becomes  small  after  only  a short  distance  across  ihe  array. 
The  actual  computer-printed  data  is  shown  in  Figures  12(b)  (l.’T  2252)  and  14(b) 

(FT  1040).  The  upper  half  of  Figure  12(b)  (FT  2252)  shows  a strong  line  at  approx- 
imately -6  Hz  Doppler  frequency  where  the  coherence  remains  relatively  high 
above  the  background  noise  level  out  to  the  maximum  distance  of  720  meters.  The 
lower  half  of  Figure  14(b)  (FT  1040)  shows  a strong  line  at  +5  Hz  with  a rapidly 
decreasing  coherence  as  the  antenna  spacing  increases.  The  coherence  is  a rela- 


tively noisy  statistical  parameter^’'  resulting  in  fluctuations  which  make  the 


8.  .lenkins,  G.  M.  , and  Watts,  D.  Cl.  (1968)  Spectral  Analysis  and  Its  Applications, 
Holden-Day,  pp  374-41  1. 

1.  Pfister,  W.  , et  al  (1975)  Pulse  Sounding  with  Closely  Spaced  Receivers  as  a 
Tool  for  Measuring  Atmospheric  Motions  and  Fine  Structure  in  the  lono- 
;,  Vol.VII,  Environmental  Research  Paper  No.  506. 
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simple  concept  of  monotonically  decreasing  coherence  with  distance  appear  to  be 
violated.  Figure  2 shows  the  least  square  linear  fit  to  the  above  illustrated  data. 
For  Flight  4-142  it  2253  UT,  an  antenna  separation  of  approximately  1450  meters 
would  be  required  for  the  coherence  to  decrease  to  a value  of  0.56,  while  for 
Flight  4-196  at  1940  UT,  the  distance  required  for  the  coherence  to  decrease  to 
0.56  is  approximately  500  meters.  As  will  be  shown  in  subsequent  reports,  this 
difference  in  behavior  can  be  related  to  the  propagation  mode;  that  is,  the  E-modes 
having  high  spatial  coherence  and  F- modes  relatively  shorter  coherence  distances. 


4.3  Aircraft  Data 

The  data  consists  of  compilation  of  samples  of  data  selected  from  six  air- 
craft missions  (Figures  3 through  14). 


5.  CONCH  SIGNS 


The  nucleus  of  this  report  comprises  compilation  of  samples  of  data  taken 
from  six  aircraft  missions  each  lasting  approximately  8 hours.  This  selection 
process  is  meant  to  illustrate  several  features  that  are  common  to  the  propagation 
of  radio  signals  from  a small  source,  basically  point-to-point  propagation  in  the 
auroral  ionosphere.  At  this  point,  it  is  necessary  to  recognize  the  fact  that  these 
aircraft  signals  are  often  spread  both  in  Doppler  frequency  and  azimuthal  arrival 
angle.  Spreading  sometimes  occurs  in  only  one  of  these  two  dimensions,  lor 
example,  spreading  occurs  in  Doppler  frequency  and  not  in  arrival  angle  as  can 
be  seen  in  Figure  14(a)  at  2038  I T.  There  are  times  when  multiple  modes  are 
set  up  which  have  different  Doppler  frequencies  and  arrival  angles  than  that  ex- 
pected from  the  known  aircraft  speed  and  direction.  A good  examp ’e  of  many 
modes  set  up  simultaneously  is  shown  in  Figure  13a  at  1719  LT  . Here  several 
modes  are  established  at  the  same  time,  very  likely  due  to  the  passage  of  a 
travelling  ionospheric  disturbance.  Another  very  common  feature  is  reflection 
from  a tilted  ionosphere  which  alters  the  arrival  angle  by  several  degrees. 

The  two  most  serious  problems  encountered  in  the  DA  ASM  experiment  at 
Goose  Bay,  in  terms  of  the  operation  of  an  over-the-horizon  radar  system  ap- 
pears to  be  TID's  and  tilted  ionospheres.  Ilow  much  of  the  occurrence  of  these 
features  can  be  ascribed  to  propagation  across  the  aurora  oval  needs  to  be  in- 
vestigated. For  this  purpose  similar  measurements  at  mid-latitude  should  be 
undertaken. 

The  relationship  between  the  observations  cataloged  here  and  ionospheric 
conditions  and  propagation  modes  will  be  presented  in  a subsequent  report. 
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Figure  12.  Flight  4-  142,  22  May  1974,  UT  2253 
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Figure  12.  Flight  4-142,  22  May  1974,  UT  0238 
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Figure  14.  Flight  4-196,  15  Jul  1974,  UT  1940 
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Figure  14.  Flight  4- 196,  15  .lul  1974,  UT  2038 
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